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COMPUTATION OF THEORETICAL BRIGHTNESS TEMPERATURES 


CORRESPONDIN'} TO THE CAPE COD CANAL RADIOMETER MEASUREMENTS 
Introduction 

The potential of microwave radiometers on eau*th satellites to sense the 
eea-surface temperature and roughness on an all-weather basis has generated in- 
terest in theoretical models for sea brightness temperatures (Tg). Wentz [l975] 
computed Tg from a two-scale scattering model and compared the results with 
Hollinger's [l97l] radiometer measurements. Although the overall agreement 
was good, discrepancies did exist and were probably due in part to errors in 
specifying the input peu*ajneters for the Tg model. Also an incomplete correction 
for the antenna patterns could have contributed to the disagreements. 

The extensive ground-truth data collected in conjunction with recent radi- 
ometer measurements of the Cape Cod Canal [Svift, 197**] and the improved de- 
termination of the permittivity of saline water [c. T. Swift, private commu- 
nication, 1975 ] now allow for a more accurate specification of the Tg model's 
input parameters. In addition. Seek [l975] has computed the corrections to Tg 
resulting from the Cape Cod Canal radiometer antenna patterns. In this report 
these investigations are used to take a closer look at the Tg model. First 
the input parameters are determined from the ground-truth data, and then the 
computations are compared with the radiometer measurements. 

It.put Parameters 

The Tg model developed by Wentz [l975] requires the following inputs. 

Operational ^drameters 

1. Radiometer frequency f 

2. Radiometer polarization P 

3 . Radiometer viewing angle 6 


I 
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Environmental Parameters 

1. Water surface temperature 

2. Water salinity a 

3. Power spectrum S(k) of the water surface, k » wavenumber 

U. Downward sky angular distribution TgCe^), 6^ ■ zenith angle 

The operational paraurters for the Cape Cod Canal measurements are f ■ l.U, 
1«.0, and 7.5 GHz; P ■ horizontal and vertical; and 6 rauiging from 20® to 70®. 
Computations are done at the three frequencies, two polarizations, and for 0 
from 0® to 70® in 10® steps . 

Two environoiental cases are considered: one between the hours of l81*0 and 

1052 EDT on Maj' 25, 1972, during which the canal's surface was rough, and the 
other between the hours of 0126 and 0130 EDT on May 26, 1972, during which the 
canal's surface was smooth. The ground-truth data required by the model are 
given in Table 1 for each case. 


Table 1. Oround-Truth Data 

Rough Surface Smooth S’orface 


Water Temperature 
Water Salinity 
Air Temperature^ 
Wind Speed ^ 


206. 1®K 
20 ®/oo 
202. 6®K 
la m/s 


203. 9®K 
20 ®/oo 
279. 8®K 
6 m/s 


if^easured approximately 58 m above the canal surface 


Permittivity of the Canal Water 

The permittivity c of saline water is a function of the water te.mpei- 
ature Ty, salinity s, and frequency f. Swift's [private communication, 1975] 
expressions for e(Ty,s,f), which are based In part on the waveguide measure- 
ments of Ho and Hall [l973] , are used, and the results are given in Table 2. 
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Table 2. Permittivity oi Canal Water 


Frequency Rough Surface 


Smooth Surface 


I.I 4 GHz 75 . UO - 50.37J 
U.O GHz 70.09 - 33.9»*J 
7.5 GHz 59.56 - 36.87J 


75.91 - >40.O2J 

69.91 - 3»*.UlJ 
58.3** - 37.73J 


Pover Spectrum of the Canal Surface 

The smooth surface was described as slick , and the wavestaff record was 
vir'.ually level. In view of this, a specular surface is assumed for the smooth- 
surface computations. The rough-surface ccaputations require that the surface 
power spectrian be specified. Information about the directional properties of 
the power spectrum of a wind-driven water surface is scarce, especially for a 
canal with an appreciable current, and the simplifying assumption that the 
spectrum is isotropic is made. Although the assumption is unrealistic, the 
anisotrophies will probably tend to average out because the Tg computations 
Involve integrations over all scatter directions. 

The isotropic power spectrum S(ic) for tc < 10"^ cm"^ was measured by a lasar 
profiloreter and by a wavestaff. However, this spectral region is unimportant 
to the Tg model because it contributes negligibly to the surface slope variance. 
The important part of the spectrum is ic > 10"^ cm“' for which no ground-truth 
spectral data is available. In order to specify S(<) in this region, the fric- 
tion velocity u« corresponding to a II 4 m/s wind in a neutral atmosphere at 58 n 
above the surface is calculated from expressions given by Cardone [1969] . The 
result, u* * 51 cm/s, is then inputted into Pierson and Stacy's [l973] expres- 
sions for the power spectrum of capillary waves (Eqs. 2.7, 2.8, and 2.9 in 
their report). In this W6iy S(ic) is determined for k 2 0.359 cm"'. One should 
note that Cardone 's expressions are for an open sea and the Pierson-Stacy ex- 
pressions are based on wave tank measurements. Applying these expressions 
to a canal possibly introduces significrnt error. 
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The capillary power spectrum is connected to the ground-truth power spec- 
trum by drawing a straight line on a log S(k) versus log k plot. This results 
in the following expressions: 


S(tc) ■ S(<j) (kj/x)** for < X < Kj (1) 

d« log[S(xj )/S(x2)]/log(x2/Xj ) (2) 

Xj - 10"3 cm"* (3) 

<2 • 0.359 cm"* (U) 

S(xj) » 1.2*10^ cm^ (laser profilometer ) (5) 

S(x2) ■ 0.372 cm^ (Pierson-Stacy ) (6) 


The spectrum for x < 10"^ cm"* can be empirically fovind from the laser profilo- 
meter data, but this is not necessary because inspection shows that this por- 
tion of the spectrum contributes less than 10"^ to the slope variar.ce. 

The treatment of two-scale scattering requires that S(x) be divided into 
a large-scale spectrum Sj^(x) and a small-scale spectrum Sg(x). 


S.(x) 



S(x) 

0 


< i IC, 


< > Kf 


Ss(<) 


0 

s(<) 



X > X 

c 


(7) 


( 8 ) 


The cutoff wavenumber is fixed by requiring that •« it ms height C of the 

small-scale spectrum equals where k denotes the radiation wavenumber. 

Ventz [1975] showed that variations in the perturbation parameter k^ from 

0.125 to 0.25 only slightly affects the 73 computations. Once x^, is found 

the slope variance of the large-scale spectrum is calculated by integrating 

x^S,(x) from x » 0 to x » x„. Table 3 gives the values of ruid 0^ for the 

c ^ c 

three frequencies. 
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Table 3« Cutoff Wavenumber and Large-r>?ale Slope Variance 
Frequency Cutoff Wavenumber »Cg Lar<?:e-Scale Slope Variance 


l.U 

GHz 

2.U<10"2 

cm"* 

1.3«10"'* 

U.O 

GHz 

3.6«10"‘ 

cm”* 

1.1*10"2 

7.5 

GHz 

7.7*10"* 

cm"* 

2.7*10"2 


Downward Sky Brightness Temperature 

The following expressions for the downward sky Tg atib-lar distribution Tg(0^) 
are used in the model: 


Ts(ez) 

(l-T)Teff ‘^'*'co8 

(9) 

T - (1 

- T /T 
•*zen'^eff- 

(10) 

Teff - 

1.12 T^r - 50 

(11) 


where is the air temperature given in Table 1, T^^^^ is the zenith atmos- 
pheric T 3 , TgQg is the cosmic Tg, and Is the air mass function [Allen, 

19 ^ 3 ] • The above equations eire the same as given by Peake [ 1969 ] except that 
the exponent kiQj.) used in place of sec 

Clear skies prevailed during the radiometer measurements, «uid Porter and 
Ventr's [l97l] values for T^gjj for a standard mid-latitude clear sky model are 
used and are given in Table U. The cosmic Tg is the sum of the galactic Tg 
and the 2.7®K resiaual Tg of the primeval fireball. The spatial variation of 
the galactic Tg is quite large. Excluding discrete sources, the galactic center 
is abcut 35 times as bright as the cold sky for the frequencies being considered 
[Allen , 1963 ] . The cold sky is assumed , and the values for the galactic Tg and 
the cosmic Tg are shown in Table h. Table 5 gives Tg(02) for the three fre- 
quencies and for the rough-surface case. The smooth-surface values are essen- 
tially the same because the 3 difference in air temperature between the two 
environmental cases results in a change in TgCo^) of only a few tenths of a K°. 
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Table U. Input Parameters for the Brightness Temperature Model 


Frequency Zenith Atmospheric Tg 


Galactic Tg Cossilc Tg 


l.U GHz 2.2®K 
U.O GHz 2.5®K 
7.5 GHz 2.9®K 


O. 50 K 3.2®K 
0.0®K 2 . 70 K 
O.QOK 2.70K 


Table 5. Downward Sky Brightness 
Zenith Angle 02 


deg 

l.U GHz 

0 

5.»* 

5 

5 . 1 * 

10 

5.U 

15 

5.5 

20 

5.5 

25 

5.6 

30 

5.T 

35 

5.9 

10 

6.0 

U 5 

6.3 

50 

6.6 

55 

T.O 

60 

7.5 

65 

8.3 

TO 

9.5 

75 

11. u 

80 

15.2 

85 

25.1 

90 

7U.U 


'emperature Angular Dist ’ . b^tion 


Tg(e2) - ®K 


U.o GHz 

7.5 GHz 

5.2 

5.6 

5.2 

5.6 

5.2 

5.6 

5.3 

5.7 

5.3 

5.8 

5.»« 

5.9 

5.6 

6.0 

5.7 

6.2 

5.9 

6.U 

6.2 

6.7 

6.5 

7.2 

7.0 

7.7 

7.6 

8.U 

8.5 

9. *4 

9.8 

11.0 

12.0 

13.5 

16.3 

18. U 

27.5 

31.2 

82.1 

92. u 


Ccrrection for Antenna Patterns 

The radiometer measurements reported by Swift [l97^] are antenna temper- 
atures T;^; that is, they are convolutions of Tg and the antenna gain patterns. 
Recent calculations by Beck [l975] show that the averaging effect of the an- 
tenna pattern produces significant differences between T;^ and Tg. These dif- 
ferences, AT » T;^ - Tg, fcr the smooth surface are given in Table 6. Calcu- 
lations were not none for 6 » 0° or IQO, and AT is assumed zero for these angles. 
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Furtheraore AT is astuned independent of the environmental parameters, and the 
values in Table 6 are also used in the rough-surface computations. The effect 
of the antenna patterns on Tg is accounted for by adding AT to the Tg model 
computations . 


Table 6. Antenna Pattern Corrections 

Horizontal Polarization 

Ta - Tg («K) 

Viewing Angle 0 


deg 

l.U GHz 

h.O GHz 

7.5 GHz 

20 

1.0 

1.0 

1.0 

30 

1.0 

1.0 

1.0 

UO 

0.5 

0.5 

1.0 

50 

1.0 

1.0 

0.5 

60 

2.0 

2.0 

1.5 

70 

3.5 

2.5 

3.5 


Vertical Polarization 


ig Angle 0 
deg 

l.U GHz 

Ta - Tg (®K) 
U.o GHz 

7.5 GHz 

20 

1.0 

0.5 

0.5 

3C 

1.5 

0.5 

1.0 

UO 

2.0 

0.5 

2.0 

50 

1.5 

0.0 

0.0 

60 

0.0 

-2.5 

-2.5 

70 

-9.0 

-8.0 

-ik.O 


results 

The model computations are compared with the measurements in Figures 1-3. 
The results for the three frequencies and for both polarizations are plotted 
versus viewing angle. The dashed and solid curves are the smooth-surface and 
rough-surface computations, respectively, and the squares and circles are the 
radiometer measurements of the smooth and rough surfaces, rci.pectively . In gen- 
eral the agreement is good especially on the change in Ta produced by the 


original 
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VIEWIN'; ANGLE 0 ( DEGREES ) 
















pretence of roughnett. Ditcreptnc^et in me^ltude do exist, but these are 
mostly iett than twice the t 2 error that is indicated by the size of the 
squares and circles. The agreement for horizontal polarization is sonewhut 
better than for vertical polarization for which the computed is greater than 
the measured at Intermediate viewing angles. The crossover of the smooth<^ 
surface and rough-surface curves for vertical polarization is in excellent agree- 
ment with tv>e ”<*asurements at 1.5 '^Hz and in fair agreement at U.O GHz. Ho 
crossover occurs for the l.U GHz computations because essentially no large- 
scale roughness is present at this f^^iquency as is indicated in Table 3. 

Possible candidates contributing to the disagreements are listed below: 

1. Error in specifying the power spectrum of the canal surface 

2. Error in specifying the permittivity of the canal water, par- 
ticularly at 1.1* GHz for which *he permittivity is a strong 
function of the water temperature and salinity 

3. • >ct of the large spatial variation of the galactic radiation 
at 1.1* GHz 

1*. Failure of the sky Tg model to consider atmospheric non-uni font- 
itles such as scattered clouds 

5. Neglect of the contribution of the canal's shore to T/^ 

6. Failure to use the total antenna pattern in computing T^ 

7. Neglect of intermediate-scale roughness by the two-scale scat- 
tering model 

0. Experimental calibration error 
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